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One of the most important questions nowadays in physics concerns the nature of the so-called 
dark energy. It is also a consensus among cosmologists that such a question will not be answered on 
the basis only of observational data. However, it is possible to diminish the range of possibilities for 
this dark component by comparing different dark energy scenarios and finding which models can be 
ruled out by current observations. In this paper, by asssuming three distinct parametrizations for 
the low-redshift evolution of the dark energy equation of state (EOS), we consider the possibility of 
discriminating between evolving dark energy and ACDM models from a joint analysis involving the 
most recent radio sources gravitational lensing sample, namely, the Cosmic All Sky Survey (CLASS) 
statistical data and the recently published gold SNe la sample. It is shown that this particular 
combination of observational data restricts considerably the dark energy parameter space, which 
enables possible distinctions between time-dependent and constant EOS's. 



PACS numbers: 98.80; 98.80.E 



I. INTRODUCTION 



The idea of a negative-pressure dark component which 
accounts for ~ 2/3 of the critical density seems to be 
strongly supported by the current cosmological observa- 
tions. Very little, however, is known about the nature of 
this extra component, a fact that has opened the possibil- 
ity for many speculations on its fundamental origin and 
has also given rise to the so-called dark energy problem. 
These speculations are usually based either on a partic- 
ular choice of the equation of state (EOS) characterizing 
the dark energy (w = p/p) or on modifications of grav- 
ity at very large scales However, irrespective of the 
choice made, a consensus is being reached nowadays on 
the fact that is extremely difficult, if not impossible, to 
single out the best-fit model of dark energy on the basis 
only of observational data (see ^ for some attempts to 
physically constrain the nature of the dark energy). 

The above conclusion, in turn, leads inevitably to a 
pessimist attitute towards the observational investigation 
of dark energy. However, as well observed by Kratochvil 
et al. Q , one may, instead of trying to find which dark 
energy scenario is correct, test which models can be ruled 
out by the available observations. This, of course, will 
not unveil the nature of the dark energy but may, with 
the current and upcoming observational data, diminish 
considerably the range of possibilitites. An interesting 
example involves two of the favorite candidates for dark 
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energy, namely, the vacuum energy density or the cosmo- 
logical constant (A) and a dynamical scalar field (0), usu- 
ally called quintessence. Among other things, what ob- 
servationally differs these two candidates for dark energy 
is that, in the former case, the EOS associated with A is 
constant along the evolution of the Universe {lo = — 1) 
whereas in generic quintessence scenarios w is a func- 
tion of the scalar field (j){t) as well as of its potential 
V{(t)) • Therefore, taking this small but important 
difference into account, one may conclude that if any 
observable deviation from a constant equation of state 
is consistently found, this naturally poses a problem for 
any model based on this assumption, which includes our 
current concordance scenario (ACDM). 

Following this reasoning, we aim, in this paper, to 
explore the prospects for constraining possible time de- 
pendence of dark energy from a joint analysis involving 
radio-selected gravitational lens statistics and supernova 
(SNe la) data. To this end, we use the most recent radio 
sources gravitational lensing sample, namely, the Cosmic 
All Sky Survey (CLASS) statistical data which consist of 
8958 radio sources out of which 13 sources are multiply 
imaged and the recently pubHshed SNe la data set 
with a total of 157 events This particular combina- 
tion of lensing statistics and SNe la data has been used 
by some authors to this very same en d la l (and also to 
place limits on constant EOS models |3, ^^^d con- 
stitutes a potential probe for possible variations of the 
dark energy EOS since it covers a considerable interval 
of redshift, which is a necessary condition to properly 
distinguish redshift-dependent equation of states from 
models with constant uj. By considering three different 
parametrizations for the time-dependence of the dark en- 
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ergy, it is shown that this combination of observational 
data restricts considerably the dark energy parameter 
space, which enables possible distinctions between time- 
dependent and constant EOS's. 



(5) 



II. uj{z) MODELS 

In this work, we are particularly interested in three 
specific parametrizations for the variation of uj with red- 
shift, i.e.. 
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where the subscript o denotes present day quantities 
and R{t) is the cosmological scale factor. The distance- 
redshift and the age-redshift relations - two fundamental 
quantities related to the observables which will be con- 
sidered in the next section - are given respectively by 



Uj{z) — UJo — UJ2\'Q-{1 + z) 



and 
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where lOq is the current value of the equation-of-state pa- 
rameter, and LOj (j = 1,2,3) are free parameters quanti- 
fying the time-dependence of the dark energy EOS, which 
must be adjusted by the observational data. Note that 
the EOS of the cosmological constant can be always re- 
covered by taking ujj = and Uo = — 1- 

The Taylor expansion (PI) was firstly suggested in Ref. 
fn\ . Constraints on (PI) were firstly studied by Cooray 
& Huterer |^ by using SNe la data, gravitational lensing 
statistics and globular clusters ages and also by Goliath 
et al. |l2l | who investigated limits to this parametrization 
from future SNe la experiments. As commented in Ref. 
0, PI is a good approximation for most quintessence 
models out to redshift of a few and it is exact for mod- 
els where the equation of state is a constant or chang- 
ing slowly. PI, however, has serious problems to explain 
age estimates of high-z objects since it predictes very 
small ages at z > 3 fl^ (In reality, PI blows up at high- 
redshifts as e^'^^^ for values of > - see Eq. (4) be- 
low). The empirical fit P2 was introduced by Efstathiou 
[T^ who argued that for a wide class of potentials asso- 
ciated to dynamical scalar field models the evolution of 
Lo{z) at z < 4 is well approximated by Eq. (2). P3 was 
recently proposed in Refs. (see also ITtII) aiming 

at solving undesirable behaviours of PI at high redshifts. 
According to Q, this parametrization is a good fit for 
many theoretically conceivable scalar field potentials, as 
well as for small recent deviations from a pure cosmolog- 
ical constant behaviour {uj = —1) [see also 0,0,123 for 
other parametrizations] . 

Since Eqs. (1-3) represent separately conserved com- 
ponents, it is straightforward to show from the energy 
conservation law [pj — —3R/R{pj + Pj)] that the ratio 
fj = Pj/ Pio foi' (P1)"(P3) evolves, respectively, as 
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where stands for the matter density parameter. 
Throughout this paper we fix r2,ji = 0.3, in accordance 
with several dynamical estimates of the quantity of mat- 
ter in the Universe |2l|. 



III. LENSING AND SNE lA CONSTRAINTS 

In our search to constrain a possible time-dependence 
of the dark energy, we adopt a joint analysis involving 
the so far largest lensing sample suitable for statistical 
analysis along with the latest SNe la data, as provided 
by Riess et al. 7]. In what follows, we discuss both 
the lensing and SNe la samples used as well as the main 
assumptions on which we performed our joint analysis. 



A. CLASS statistical Sample 

The final CLASS well-defined statistical sample con- 
sists of 8958 radio sources out of which 13 sources are 
multiply imaged. Here we work only with those mul- 
tiply imaged sources whose image-splittings are known 
(or likely) to be caused by single galaxies. There are 
9 such radio sources: 0218-f357, 0445-^123, 0631-1-519, 
0712-1-472, 0850-^054, 1152-1-199, 1422-1-231, 1933-1-503, 
2319-f051. We thus work with a total of 8954 radio 
sources "2^, 0- The sources probed by CLASS at 
v = 5 GHz are well represented by power-law differential 
number-flux density relation: \dN/dS\ cx {S/S'^)^ with 
7] = 2.07 ± 0.02 (1.97 ± 0.14) for S > S" {< S°) where 
— 30 mJy |2J|. The CLASS unlensed sources can 
be adequately described by a Gaussian model with mean 
redshift, z = 1.27 and a dispersion of 0.95. 

We start our analysis by assuming the singular isother- 
mal sphere (SIS) model for the lens mass distribution. As 
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has been discussed elsewhere this assumption represents 
a good approximation to the real mass distribution in 
galaxies (see, e.g., H^). For the present analysis we also 
ignore the evolution of the number density of galaxies and 
assume that the comoving number density is conserved. 
The present day comoving number density of galaxies is 



(j){L)dL, 



(9) 



where^(L) is the well known Schechter Luminosity Func- 
tion |2fl |. 

The differential optical depth of lensing in traversing 
dzL with angular separation between ip and ip+dcp is ,2 7.] : 
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where the function F* is defined as 



p* 
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The quantities Dql, Dqs and Dls represent, respec- 
tively, the angular diameter distances from the observer 
to the lens, from the observer to the source and between 
the lens and the source. In order to relate the characteris- 
tic luminosity L* to the characteristic velocity dispersion 
, we use the Faber- Jackson relation for early- type 
galaxies (L* oc u*'''), with 7 = 4. For the analysis pre- 
sented here we neglect the contribution of spirals as lenses 
because their velocity dispersion is small when compared 
to ellipticals. 

The two large-scale galaxy surveys, namely, the 2dF- 
GRS (3^] and the SDSS 35] have produced converg- 
ing results on the total LF. The surveys determined the 
Schechter parameters for galaxies (all types) at z < 0.2. 
Chae (2^] has worked extensively on the information pro- 
vided by these recent galaxy surveys to extract the lo- 
cal type-specific LFs. For our analysis here, we adopt 
the normalization corrected Schechter parameters of the 
2dFGRS survey jH l^: a = -0.74, 0* = 0.82 x 
IQ-'^h^Mpc-^, V* = 185km/s and F* = 0.014. 

The normalized image angular separation distribution 
for a source at 2:5 is obtained by integrating ^^^'^^ over 

ZL- 
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The corrected (for magnification and selection effects) 
image separation distribution function for a single source 



at redshift zs is given by [3C 
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Similarly, the corrected lensing probability for a given 
source at redshift z is given by 



(14) 



Here (h and AO are related to as d — „ and B 

is the magnification bias. This is considered because, as 
is widely known, gravitational lensing causes a magni- 
fication of images and this transfers the lensed sources 
to higher flux density bins. In other words, the lensed 
sources are over-represented in a flux-limited sample. 
The magnification bias B{zs, Su) increases the lensing 
probability significantly in a bin of total fiux density {S^) 
by a factor 



dN,,i>S,) 



dSy 



(15) 
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In the above expression, Nzs{> S„) is the intrinsic flux 
density relation for the source population at redshift zs- 
Nzs (> S^) gives the number of sources at redshift zs hav- 
ing flux greater than 5^. For the SIS model, the mag- 
nification probability distribution is p(/x) = 8/^'^. The 
minimum and maximum total magnifications ^mm and 
Mniaa: in equation Hlt)|l depend on the observational char- 
acteristics as well as on the lens model. For the SIS 
model, the minimum total magnification is ^min — 2 and 
the maximum total magnification is /i,noa; = 00. The 
magnification bias B depends on the differential number- 
flux density relation \dNzs{> S^)/dS,y\. The differential 
number-flux relation needs to be known as a function of 
the source redshift. At present redshifts of only a few 
CLASS sources are known. We, therefore, ignore red- 
shift dependence of the differential number-flux density 
relation. Following Chae |22|, we further ignore the de- 
pendence of the differential number-flux density relation 
on the spectral index of the source. 

Two important selection criteria for CLASS statisti- 
cal sample are (i) that the ratio of the flux densities of 
the fainter to the brighter images TZmin is > 0.1. Given 
such an observational limit, the minimum total magnifi- 
cation for double imaging the adopted model of the lens 
is f^min = 2(1 4- TZ„iin/l - TZmin) [111; (ii) that the im- 
age components in lens systems must have separations 
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FIG. 1: Current constraints on the plane u — ujj from CLASS lensing statistics and SNe la data. Note that for the three 
different parametrizations considered in this paper the joint CLASS + SNe la analysis clearly prefers the regions where ujj 7^ 0, 
i.e., a time-dependent EOS. a) Confidence regions (68.3% and 95.4%) in the plane u — u\ for PI. b) The same as in Panel la 
for P2. c) The same as in Panel la for P3. 



> 0.3 arcsec. We incorporate this selection criterion by 
setting the lower limit of A6' in equation H14|) as 0.3 arc- 
sec. 

B. SNe la sample. 

The SNe la sample of Riess et al. consists of 186 
events distributed over the redshift interval 0.01 < z < 
1.7 and constitutes the compilation of the best observa- 
tions made so far by the two supernova search teams plus 
16 new events observed by HST. This total data-set was 
divided into "high-confidence" {gold) and "likely but not 
certain" [silver) subsets. Here, we will consider only the 
157 events that constitute the so-called gold sample. The 
best fit to the set of parameters s is obtained by using a 
statistics, i.e., 

2 [Mp(2:|s) -Mo(2:|s)]^ 

XsNe = 2^ -2 ' (16) 

1=1 i 

where iJLp{z\s) = m — M = Slogd^ + 25 is the predicted 
distance modulus for a supernova at redshift z, /i*(z|s) 
is the extinction corrected distance modulus for a given 
SNe la at z^, and Ui is the uncertainty in the individual 
distance moduli, which includes uncertainties in galaxy 
redshift due to a peculiar velocity of 400 km/s. The Hub- 
ble parameter Ho is considered a "nuisance" parameter 
so that we marginalize over it (For some recent SNe la 
studies, see jsj)- 



C. Results. 

The main results of our joint analysis are displayed 
in Figure 1. Panels la-lc show the confidence regions 
(68.3% and 95.4%) in the ujo - ojj plane for PI, P2 and 
P3, respectively. As stated earlier, the matter density 
parameter has been fixed in all analyses at = 0.3, 
in agreement with some clustering estimates l2lj. The 
contours are defined by the conventional two-parameter 
levels (2.30 and 6.17), where Xjoj^/ — XsNe ~ ^l^i^ ^^id 

1 = Ciens/^l^x is the normahzed hkelihood for lenses. 
From this combination of observational data we find 
that the best-fit parameters for PI are ujo = —1.5 and 
LUi = 2.1, which corresponds to an accelerating scenario 
with transition redshift (at which the Universe switches 
from acceleration to acceleration) zt ~ 0.26 and a total 
expanding age of ~ 7.0h~^ Gyr. For parametrizations 

2 and 3 the best-fit scenarios occur at ojo = —1.4 and 
UJ2 — —2 and Uq = —1.7 and uj^ = 4.1, corresponding to 
a 10/i~^-Gyr-old and 7.5ft.^^-Gyr-old universes with tran- 
sition redshifts at Zt — 0.53 and Zt — 0.31, respectively. 
Note that the estimates of Zt for P2 and P3 are inside 
the 1(7 interval inferred for the transition redshift from 
the current SNe la data (see also [l^). The best-fit 
values for each parametrization also imply a beginning of 
a phantom behavior at Zph — 0.23 (PI), Zph — 0.22 (P2) 
and Zph ~ 0.20 (P3). 

From the above figures, however, it is clear the most 
important conclusion one may reach resides on the fact 
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that for the three parametrizations considered in this pa- 
per the joint CLASS + SNe la analysis clearly prefers 
the regions where luj ^ 0, i.e., a time-dependent EOS. 
In particular, note that the point cuj = 0, which means a 
time-independent EOS, is at least 2a off from the central 
values obtained for this parameter. Such a result seems 
to be more restrictive than that obtained in Ref. (sse 
also "20]), in which a combination of the latest SNe la, 
cosmic microwave background and large-scale structure 
data showed no hint of departures from the model corre- 
sponding to Einstein's cosmological constant (ujj = and 
Wo = — 1). It is possible that this particular preference 
for a time-dependent EOS changes if the current cosmic 
microwave background (CMB) data were included in the 
analysis. As discussed in |23, an evolving dark energy 
EOS affects the features of temperature anisotropics in 
CMB in at least two ways, namely, the position of the 
acoustic peaks as well as the integrated Sachs-Wolf efect. 
However, if independent analyses involving a more signif- 
icant number of different data sets confirm this preference 
of the observational data for values of uij 7^ 0, these re- 
sults surely will bring to light a new consistency problem 
for our current standard cosmological model since for this 
case ujj is necessarily null. 



of high-quality observational data. Surely, the most re- 
markable among these concepts is the idea of a dark 
energy-dominated universe, which is motivated from an 
impressive convergence of independent observational re- 
sults. This idea in turn gave rise to the so-called dark 
energy problem since the nature of this dark component 
is completely unknown at present. 

In this paper, although aware of the impossibility of 
determining the nature of the dark energy on the basis 
only of observational data, we have considered the pos- 
sibility of discriminating two of our favorite candidates 
for this mysterious component, namely, the cosmological 
constant (A) and a dynamical scalar field (</)) [33]. By 
considering three different parametrizations for the dark 
enegy EOS, we have placed limits on the time-dependent 
term of these parametrizations (wj) from a joint analy- 
sis involving the most recent radio sources gravitational 
lensing sample and SNe la data. We have shown that 
this particular combination of observational data prefer 
values of Wj ^ 0, i.e., a time-dependent EOS. We believe 
that if such a result is confirmed by the upcoming obser- 
vations, it may shed some light on our search for a better 
understanding of the nature of the so-called dark energy. 
Acknowledgments 



IV. CONCLUSION 

Very recently, the field of cosmology has entered a 
golden era. An era where new and revolutionary con- 
cepts were introduced with the support of a plethora 
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